An integrated technology for quantitative wide mutational scanning of human antibody Fab libraries.

Antibodies are engineerable quantities in medicine. Learning antibody molecular recognition would enable the in silico design of high affinity binders against nearly any proteinaceous surface. Yet, publicly available experiment antibody sequence-binding datasets may not contain the mutagenic, antigenic, or antibody sequence diversity necessary for deep learning approaches to capture molecular recognition. In part, this is because limited experimental platforms exist for assessing quantitative and simultaneous sequence-function relationships for multiple antibodies. Here we present MAGMA-seq, an integrated technology that combines multiple antigens and multiple antibodies and determines quantitative biophysical parameters using deep sequencing. We demonstrate MAGMA-seq on two pooled libraries comprising mutants of nine different human antibodies spanning light chain gene usage, CDR H3 length, and antigenic targets. We demonstrate the comprehensive mapping of potential antibody development pathways, sequence-binding relationships for multiple antibodies simultaneously, and identification of paratope sequence determinants for binding recognition for broadly neutralizing antibodies (bnAbs). MAGMA-seq enables rapid and scalable antibody engineering of multiple lead candidates because it can measure binding for mutants of many given parental antibodies in a single experiment.

Computational Modeling and Characterization of Peptides Derived from Nanobody Complementary-Determining Region 2 (CDR2) Targeting Active-State Conformation of the ß2-Adrenergic Receptor (ß2AR).

This study assessed the suitability of the complementarity-determining region 2 (CDR2) of the nanobody (Nb) as a template for the derivation of nanobody-derived peptides (NDPs) targeting active-state ß2-adrenergic receptor (ß2AR) conformation. Sequences of conformationally selective Nbs favoring the agonist-occupied ß2AR were initially analyzed by the informational spectrum method (ISM). The derived NDPs in complex with ß2AR were subjected to protein-peptide docking, molecular dynamics (MD) simulations, and metadynamics-based free-energy binding calculations. Computational analyses identified a 25-amino-acid-long CDR2-NDP of Nb71, designated P4, which exhibited the following binding free-energy for the formation of the ß2AR:P4 complex (ΔG = -6.8 ± 0.8 kcal/mol or a Ki = 16.5 µM at 310 K) and mapped the ß2AR:P4 amino acid interaction network. In vitro characterization showed that P4 (i) can cross the plasma membrane, (ii) reduces the maximum isoproterenol-induced cAMP level by approximately 40% and the isoproterenol potency by up to 20-fold at micromolar concentration, (iii) has a very low affinity to interact with unstimulated ß2AR in the cAMP assay, and (iv) cannot reduce the efficacy and potency of the isoproterenol-mediated ß2AR/ß-arrestin-2 interaction in the BRET2-based recruitment assay. In summary, the CDR2-NDP, P4, binds preferentially to agonist-activated ß2AR and disrupts Gαs-mediated signaling.

Do antibody CDR loops change conformation upon binding?

Antibodies have increasingly been developed as drugs with over 100 now licensed in the US or EU. During development, it is often necessary to increase or reduce the affinity of an antibody and rational attempts to do so rely on having a structure of the antibody-antigen complex often obtained by modeling. The antigen-binding site consists primarily of six loops known as complementarity-determining regions (CDRs), and an open question has been whether these loops change their conformation when they bind to an antigen. Existing surveys of antibody-antigen complex structures have only examined CDR conformational change in case studies or small-scale surveys. With an increasing number of antibodies where both free and complexed structures have been deposited in the Protein Data Bank, a large-scale survey of CDR conformational change during binding is now possible. To this end, we built a dataset, AbAgDb, that currently includes 177 antibodies with high-quality CDRs, each of which has at least one bound and one unbound structure. We analyzed the conformational change of the Cα backbone of each CDR upon binding and found that, in most cases, the CDRs (other than CDR-H3) show minimal movement, while 70.6% and 87% of CDR-H3s showed global Cα RMSD ≤ 1.0Å and ≤ 2.0Å, respectively. We also compared bound CDR conformations with the conformational space of unbound CDRs and found most of the bound conformations are included in the unbound conformational space. In future, our results will contribute to developing insights into antibodies and new methods for modeling and docking.

AbFlex: designing antibody complementarity determining regions with flexible CDR definition.

MOTIVATION: Antibodies are proteins that the immune system produces in response to foreign pathogens. Designing antibodies that specifically bind to antigens is a key step in developing antibody therapeutics. The complementarity determining regions (CDRs) of the antibody are mainly responsible for binding to the target antigen, and therefore must be designed to recognize the antigen. RESULTS: We develop an antibody design model, AbFlex, that exhibits state-of-the-art performance in terms of structure prediction accuracy and amino acid recovery rate. Furthermore, >38% of newly designed antibody models are estimated to have better binding energies for their antigens than wild types. The effectiveness of the model is attributed to two different strategies that are developed to overcome the difficulty associated with the scarcity of antibody-antigen complex structure data. One strategy is to use an equivariant graph neural network model that is more data-efficient. More importantly, a new data augmentation strategy based on the flexible definition of CDRs significantly increases the performance of the CDR prediction model. AVAILABILITY AND IMPLEMENTATION: The source code and implementation are available at https://github.com/wsjeon92/AbFlex.

Investigating the ability of deep learning-based structure prediction to extrapolate and/or enrich the set of antibody CDR canonical forms.

Deep learning models have been shown to accurately predict protein structure from sequence, allowing researchers to explore protein space from the structural viewpoint. In this paper we explore whether "novel" features, such as distinct loop conformations can arise from these predictions despite not being present in the training data. Here we have used ABodyBuilder2, a deep learning antibody structure predictor, to predict the structures of ~1.5M paired antibody sequences. We examined the predicted structures of the canonical CDR loops and found that most of these predictions fall into the already described CDR canonical form structural space. We also found a small number of "new" canonical clusters composed of heterogeneous sequences united by a common sequence motif and loop conformation. Analysis of these novel clusters showed their origins to be either shapes seen in the training data at very low frequency or shapes seen at high frequency but at a shorter sequence length. To evaluate explicitly the ability of ABodyBuilder2 to extrapolate, we retrained several models whilst withholding all antibody structures of a specific CDR loop length or canonical form. These "starved" models showed evidence of generalisation across CDRs of different lengths, but they did not extrapolate to loop conformations which were highly distinct from those present in the training data. However, the models were able to accurately predict a canonical form even if only a very small number of examples of that shape were in the training data. Our results suggest that deep learning protein structure prediction methods are unable to make completely out-of-domain predictions for CDR loops. However, in our analysis we also found that even minimal amounts of data of a structural shape allow the method to recover its original predictive abilities. We have made the ~1.5 M predicted structures used in this study available to download at https://doi.org/10.5281/zenodo.10280181.

Structure- and machine learning-guided engineering demonstrate that a non-canonical disulfide in an anti-PD-1 rabbit antibody does not impede antibody developability.

Rabbits produce robust antibody responses and have unique features in their antibody repertoire that make them an attractive alternative to rodents for in vivo discovery. However, the frequent occurrence of a non-canonical disulfide bond between complementarity-determining region (CDR) H1 (C35a) and CDRH2 (C50) is often seen as a liability for therapeutic antibody development, despite limited reports of its effect on antibody binding, function, and stability. Here, we describe the discovery and humanization of a human-mouse cross-reactive anti-programmed cell death 1 (PD-1) monoclonal rabbit antibody, termed h1340.CC, which possesses this non-canonical disulfide bond. Initial removal of the non-canonical disulfide resulted in a loss of PD-1 affinity and cross-reactivity, which led us to explore protein engineering approaches to recover these. First, guided by the sequence of a related clone and the crystal structure of h1340.CC in complex with PD-1, we generated variant h1340.SA.LV with a potency and cross-reactivity similar to h1340.CC, but only partially recovered affinity. Side-by-side developability assessment of both h1340.CC and h1340.SA.LV indicate that they possess similar, favorable properties. Next, and prompted by recent developments in machine learning (ML)-guided protein engineering, we used an unbiased ML- and structure-guided approach to rapidly and efficiently generate a different variant with recovered affinity. Our case study thus indicates that, while the non-canonical inter-CDR disulfide bond found in rabbit antibodies does not necessarily constitute an obstacle to therapeutic antibody development, combining structure- and ML-guided approaches can provide a fast and efficient way to improve antibody properties and remove potential liabilities.

Matrixed CDR grafting: A neoclassical framework for antibody humanization and developability.

Discovery and optimization of a biotherapeutic monoclonal antibody requires a careful balance of target engagement and physicochemical developability properties. To take full advantage of the sequence diversity provided by different antibody discovery platforms, a rapid and reliable process for humanization of antibodies from nonhuman sources is required. Canonically, maximizing homology of the human variable region (V-region) to the original germline was believed to result in preservation of binding, often without much consideration for inherent molecular properties. We expand on this approach by grafting the complementary determining regions (CDRs) of a mouse anti-LAG3 antibody into an extensive matrix of human variable heavy chain (VH) and variable light chain (VL) framework regions with substantially broader sequence homology to assess the impact on complementary determining region-framework compatibility through progressive evaluation of expression, affinity, biophysical developability, and function. Specific VH and VL framework sequences were associated with major expression and purification phenotypes. Greater VL sequence conservation was correlated with retained or improved affinity. Analysis of grafts that bound the target demonstrated that initial developability criteria were significantly impacted by VH, but not VL. In contrast, cell binding and functional characteristics were significantly impacted by VL, but not VH. Principal component analysis of all factors identified multiple grafts that exhibited more favorable antibody properties, notably with nonoptimal sequence conservation. Overall, this study demonstrates that modern throughput systems enable a more thorough, customizable, and systematic analysis of graft-framework combinations, resulting in humanized antibodies with improved global properties that may progress through development more quickly and with a greater probability of success.

IGL CDR3 Hydropathy and Antigen Chemical Complementarity Associated with Greater Disease-Free Survival in Lung Adenocarcinoma: Implications for Gender Disparities.

With lung cancer remaining a challenging disease, new approaches to biomarker discovery and therapy development are needed. Recent immunogenomics, adaptive immune receptor approaches have indicated that it is very likely that B cells play an important role in mediating better overall outcomes. As such, we assessed physicochemical features of lung adenocarcinoma resident IGL complementarity determining region-3 (CDR3) amino acid (AA) sequences and determined that hydrophobic CDR3 AA sequences were associated with a better disease-free survival (DFS) probability. Further, using a recently developed chemical complementarity scoring algorithm particularly suitable for the evaluation of large patient datasets, we determined that IGL CDR3 chemical complementarity with certain cancer testis antigens was associated with better DFS. Chemical complementarity scores for IGL CDR3-MAGEC1 represented a gender bias, with an overrepresentation of males among the higher IGL-CDR3-CTA complementarity scores that were in turn associated with better DFS (logrank p < 0.065). Overall, this study pointed towards potential biomarkers for prognoses that, in some cases are likely gender-specific; and towards biomarkers for guiding therapy, e.g., IGL-based opportunities for antigen targeting in the lung cancer setting.

Computational optimization of antibody humanness and stability by systematic energy-based ranking.

Conventional methods for humanizing animal-derived antibodies involve grafting their complementarity-determining regions onto homologous human framework regions. However, this process can substantially lower antibody stability and antigen-binding affinity, and requires iterative mutational fine-tuning to recover the original antibody properties. Here we report a computational method for the systematic grafting of animal complementarity-determining regions onto thousands of human frameworks. The method, which we named CUMAb (for computational human antibody design; available at http://CUMAb.weizmann.ac.il ), starts from an experimental or model antibody structure and uses Rosetta atomistic simulations to select designs by energy and structural integrity. CUMAb-designed humanized versions of five antibodies exhibited similar affinities to those of the parental animal antibodies, with some designs showing marked improvement in stability. We also show that (1) non-homologous frameworks are often preferred to highest-homology frameworks, and (2) several CUMAb designs that differ by dozens of mutations and that use different human frameworks are functionally equivalent.

abYpap: improvements to the prediction of antibody VH/VL packing using gradient boosted regression.

The Fv region of the antibody (comprising VH and VL domains) is the area responsible for target binding and thus the antibody's specificity. The orientation, or packing, of these two domains relative to each other influences the topography of the Fv region, and therefore can influence the antibody's binding affinity. We present abYpap, an improved method for predicting the packing angle between the VH and VL domains. With the large data set now available, we were able to expand greatly the number of features that could be used compared with our previous work. The machine-learning model was tuned for improved performance using 37 selected residues (previously 13) and also by including the lengths of the most variable 'complementarity determining regions' (CDR-L1, CDR-L2 and CDR-H3). Our method shows large improvements from the previous version, and also against other modeling approaches, when predicting the packing angle.

Conformational features and interaction mechanisms of VH H antibodies with ß-hairpin CDR3: A case of Nb8-HigB2 interaction.

The ß-hairpin conformation is regarded as an important basic motif to form and regulate protein-protein interactions. Single-domain VH H antibodies are potential therapeutic and diagnostic tools, and the third complementarity-determining regions of the heavy chains (CDR3s) of these antibodies are critical for antigen recognition. Although the sequences and conformations of the CDR3s are diverse, CDR3s sometimes adopt ß-hairpin conformations. However, characteristic features and interaction mechanisms of ß-hairpin CDR3s remain to be fully elucidated. In this study, we investigated the molecular recognition of the anti-HigB2 VH H antibody Nb8, which has a CDR3 that forms a ß-hairpin conformation. The interaction was analyzed by evaluation of alanine-scanning mutants, molecular dynamics simulations, and hydrogen/deuterium exchange mass spectrometry. These experiments demonstrated that positions 93 and 94 (Chothia numbering) in framework region 3, which is right outside CDR3 by definition, play pivotal roles in maintaining structural stability and binding properties of Nb8. These findings will facilitate the design and optimization of single-domain antibodies.

LC-MS Approach to Decipher a Light Chain Chromatographic Peak Splitting of a Monoclonal Antibody.

PURPOSE: Monoclonal antibodies (mAbs), like other protein therapeutics, are prone to various forms of degradation, some of which are difficult to distinguish from the native form yet may alter potency. A generalizable LC-MS approach was developed to enable quantitative analysis of isoAsp. In-depth understanding of product quality attributes (PQAs) enables optimization of the manufacturing process, better formulation selection, and decreases risk associated with product handling in the clinic or during shipment. METHODS: Reversed-phase chromatographic peak splitting was observed when a mAb was exposed to elevated temperatures. Multiple LC-MS based methods were applied to identify the reason for peak splitting. The approach involved the use of complementary HPLC columns, multiple enzymatic digestions and different MS/MS ion dissociation methods. In addition, mAb potency was measured by enzyme-linked immunosorbent assay (ELISA). RESULTS: The split peaks had identical masses, and the root cause of the peak splitting was identified as isomerization of an aspartic acid located in the complementarity-determining region (CDR) of the light chain. And the early eluting and late eluting peaks were collected and performed enzymatic digestion to confirm the isoAsp enrichment in the early eluting peak. In addition, decreased potency was observed in the same heat-stressed sample, and the increased isoAsp levels in the CDR correlate well with a decrease of potency. CONCLUSION: Liquid chromatography-mass spectrometry (LC-MS) has been utilized extensively to assess PQAs of biological therapeutics. In this study, a generalizable LC-MS-based approach was developed to enable identification and quantitation of the isoAsp-containing peptides.

A window into the human immune system: comprehensive characterization of the complexity of antibody complementary-determining regions in functional antibodies.

The human immune system uses antibodies to neutralize foreign antigens. They are composed of heavy and light chains, both with constant and variable regions. The variable region has six hypervariable loops, also known as complementary-determining regions (CDRs) that determine antibody diversity and antigen specificity. Knowledge of their significance, and certain residues present in these areas, is vital for antibody therapeutics development. This study includes an analysis of more than 11,000 human antibody sequences from the International Immunogenetics information system (IMGT). The analysis included parameters such as length distribution, overall amino acid diversity, amino acid frequency per CDR and residue position within antibody chains. Overall, our findings confirm existing knowledge, such as CDRH3's high length diversity and amino acid variability, increased aromatic residue usage, particularly tyrosine, charged and polar residues like aspartic acid, serine, and the flexible residue glycine. Specific residue positions within each CDR influence these occurrences, implying a unique amino acid type distribution pattern. We compared amino acid type usage in CDRs and non-CDR regions, both in globular and transmembrane proteins, which revealed distinguishing features, such as increased frequency of tyrosine, serine, aspartic acid, and arginine. These findings should prove useful for future optimization, improvement of affinity, synthetic antibody library design, or the creation of antibodies de-novo in silico.

Benchmarking TriadAb using targets from the second antibody modeling assessment.

Computational modeling and design of antibodies has become an integral part of today's research and development in antibody therapeutics. Here we describe the Triad Antibody Homology Modeling (TriadAb) package, a functionality of the Triad protein design platform that predicts the structure of any heavy and light chain sequences of an antibody Fv domain using template-based modeling. To gauge the performance of TriadAb, we benchmarked against the results of the Second Antibody Modeling Assessment (AMA-II). On average, TriadAb produced main-chain carbonyl root-mean-square deviations between models and experimentally determined structures at 1.10 Å, 1.45 Å, 1.41 Å, 3.04 Å, 1.47 Å, 1.27 Å, 1.63 Å in the framework and the six complementarity-determining regions (H1, H2, H3, L1, L2, L3), respectively. The inaugural results are comparable to those reported in AMA-II, corroborating with our internal bench-based experiences that models generated using TriadAb are sufficiently accurate and useful for antibody engineering using the sequence design capabilities provided by Triad.

A library approach for the de novo high-throughput isolation of humanized VHH domains with favorable developability properties following camelid immunization.

In this study, we generated a novel library approach for high throughput de novo identification of humanized single-domain antibodies following camelid immunization. To achieve this, VHH-derived complementarity-determining regions-3 (CDR3s) obtained from an immunized llama (Lama glama) were grafted onto humanized VHH backbones comprising moderately sequence-diversified CDR1 and CDR2 regions similar to natural immunized and naïve antibody repertoires. Importantly, these CDRs were tailored toward favorable in silico developability properties, by considering human-likeness as well as excluding potential sequence liabilities and predicted immunogenic motifs. Target-specific humanized single-domain antibodies (sdAbs) were readily obtained by yeast surface display. We demonstrate that, by exploiting this approach, high affinity sdAbs with an optimized in silico developability profile can be generated. These sdAbs display favorable biophysical, biochemical, and functional attributes and do not require any further sequence optimization. This approach is generally applicable to any antigen upon camelid immunization and has the potential to significantly accelerate candidate selection and reduce risks and attrition rates in sdAb development.

TRB CDR3-cancer testis antigen chemical complementarity scoring for identifying productive immune responses in renal cell carcinoma.

BACKGROUND: Immunogenomics approaches to the characterization of renal cell carcinoma (RCC) have helped to better our understanding of the features of RCC immune dysfunction. However, much is still unknown with regard to specific immune interactions and their impact in the tumor microenvironment. OBJECTIVE: This study applied chemical complementarity scoring for the TRB complementarity determining region-3 (CDR3) amino acid sequences and cancer testis antigens (CTAs) to determine whether such complementarity correlated with survival and the expression of immune marker genes. METHODS: TRB recombination reads from RCC tumor samples from RNAseq files obtained from two separate databases, Moffitt Cancer Center and The Cancer Genome Atlas (TCGA), were evaluated. Chemical complementarity scores (CSs) were calculated for TRB CDR3-CTA pairs and survival assessments based on those CSs were performed. RESULTS: Moffitt Cancer Center and TCGA cases representing the upper 50th percentile of chemical CSs for TRB CDR3 amino acid sequences and the CTA POTEA were found to be associated with a better overall survival (OS) Also, greater tumor RNA expression of multiple immune signature genes, including granzyme A, granzyme B, and interferon-gamma were correlated with the higher chemical CSs. CONCLUSIONS: These results indicate that TRB CDR3-CTA chemical complementarity scoring may be useful in distinguishing RCC cases with a productive, anti-tumor immune response from cases where basic immune parameter assessments are inconsistent with a productive immune response.

Arginine cluster introduction on framework region in anti-lysozyme antibody improved association rate constant by changing conformational diversity of CDR loops.

Antibodies are used for many therapeutic and biotechnological purposes. Because the affinity of an antibody to the antigen is critical for clinical efficacy of pharmaceuticals, many affinity maturation strategies have been developed. Although we previously reported an affinity maturation strategy in which the association rate of the antibody toward its antigen is improved by introducing a cluster of arginine residues into the framework region of the antibody, the detailed molecular mechanism responsible for this improvement has been unknown. In this study, we introduced five arginine residues into an anti-hen egg white lysozyme antibody (HyHEL10) Fab fragment to create the R5-mutant and comprehensively characterized the interaction between antibody and antigen using thermodynamic analysis, X-ray crystallography, and molecular dynamics (MD) simulations. Our results indicate that introduction of charged residues strongly enhanced the association rate, as previously reported, and the antibody-antigen complex structure was almost the same for the R5-mutant and wild-type Fabs. The MD simulations indicate that the mutation increased conformational diversity in complementarity-determining region loops and thereby enhanced the association rate. These observations provide the molecular basis of affinity maturation by R5 mutation.

Blood-based T cell receptor anti-viral CDR3s are associated with worse overall survival for neuroblastoma.

With the advent of large collections of adaptive immune receptor recombination reads representing cancer, there is the opportunity to further investigate the adaptive immune response to viruses in the cancer setting. This is a particularly important goal due to longstanding but still not well-resolved questions about viral etiologies in cancer and viral infections as comorbidities. In this report, we assessed the T cell receptor complementarity determining region-3 (CDR3) amino acid (AA) sequences, for blood-sourced TCRs from neuroblastoma (NBL) cases, for exact AA sequence matches to previously identified anti-viral TCR CDR3 AA sequences. Results indicated the presence of anti-viral TCR CDR3 AA sequences in the NBL blood samples highly significantly correlated with worse overall survival. Furthermore, the TCR CDR3 AA sequences demonstrating chemical complementarity to many cytomegalovirus antigens represented cases with a worse outcome, including cases where such CDR3s were obtained from tumor samples. Overall, these results indicate a significant need for, and provide a novel strategy for assessing viral infection complications in NBL patients.

KA-Search, a method for rapid and exhaustive sequence identity search of known antibodies.

Antibodies with similar amino acid sequences, especially across their complementarity-determining regions, often share properties. Finding that an antibody of interest has a similar sequence to naturally expressed antibodies in healthy or diseased repertoires is a powerful approach for the prediction of antibody properties, such as immunogenicity or antigen specificity. However, as the number of available antibody sequences is now in the billions and continuing to grow, repertoire mining for similar sequences has become increasingly computationally expensive. Existing approaches are limited by either being low-throughput, non-exhaustive, not antibody specific, or only searching against entire chain sequences. Therefore, there is a need for a specialized tool, optimized for a rapid and exhaustive search of any antibody region against all known antibodies, to better utilize the full breadth of available repertoire sequences. We introduce Known Antibody Search (KA-Search), a tool that allows for the rapid search of billions of antibody variable domains by amino acid sequence identity across either the variable domain, the complementarity-determining regions, or a user defined antibody region. We show KA-Search in operation on the [Formula: see text]2.4 billion antibody sequences available in the OAS database. KA-Search can be used to find the most similar sequences from OAS within 30 minutes and a representative subset of 10 million sequences in less than 9 seconds. We give examples of how KA-Search can be used to obtain new insights about an antibody of interest. KA-Search is freely available at https://github.com/oxpig/kasearch .

Understanding the complementarity and plasticity of antibody-antigen interfaces.

MOTIVATION: While antibodies have been ground-breaking therapeutic agents, the structural determinants for antibody binding specificity remain to be fully elucidated, which is compounded by the virtually unlimited repertoire of antigens they can recognize. Here, we have explored the structural landscapes of antibody-antigen interfaces to identify the structural determinants driving target recognition by assessing concavity and interatomic interactions. RESULTS: We found that complementarity-determining regions utilized deeper concavity with their longer H3 loops, especially H3 loops of nanobody showing the deepest use of concavity. Of all amino acid residues found in complementarity-determining regions, tryptophan used deeper concavity, especially in nanobodies, making it suitable for leveraging concave antigen surfaces. Similarly, antigens utilized arginine to bind to deeper pockets of the antibody surface. Our findings fill a gap in knowledge about the antibody specificity, binding affinity, and the nature of antibody-antigen interface features, which will lead to a better understanding of how antibodies can be more effective to target druggable sites on antigen surfaces. AVAILABILITY AND IMPLEMENTATION: The data and scripts are available at: https://github.com/YoochanMyung/scripts.